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Abstract 



A constituent parton picture of hadrons with logarithmic confinement arises natu- 
rally in light-front QCD when the hamiltonian is computed using a perturbative renor- 
malization group; however, rotational symmetry is not manifest as a simple kinematic 
symmetry. Relevant operators must typically be fine-tuned when a perturbative renor- 
malization group is used, and in light-front field theory such operators contain functions 
of longitudinal momenta. We explore the possibility that a gluon mass operator with sim- 
ple longitudinal momentum dependence may improve approximate kinematic rotational 
symmetry, as revealed by the potential between heavy quarks. 
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The solution of Quantum Chromodynamics in the nonperturbative domain remains 
one of the most important problems in physics. A nonperturbative approach based on 
light-front quantization has been suggested by Wilson et al Q. The basic assumption is 
that it is possible to derive a constituent picture for hadrons from QCD [||, ^, ^. If this is 
possible, nonperturbative bound state problems in QCD can be approximated as coupled, 
few-body Schrodinger equations. 

To arrive at a constituent approximation, vacuum fluctuations need to be isolated. 
This is achieved by formulating the theory in light-front coordinates, in which the vacuum 
is made trivial by the introduction of a longitudinal momentum cutoff and vacuum effects 
must be replaced by effective interactions. However, a trivial vacuum alone does not 
ensure the emergence of a constituent picture. There may still be infinitely many Fock 
components needed to describe a hadron. In addition to removing vacuum components 
from the states, significant coupling between few-body and many-body states must be 
replaced by effective interactions involving few-body states alone. This can be achieved, 
if the many-body states have high energy, by a renormalization group Q . 

We sketch the simplest version of this procedure and motivate the current calculations. 
The first step is to use a similarity renormalization group Q and coupling coherence Q 
to compute the renormalized hamiltonian to a fixed order in the running canonical cou- 
pling. The truncation of this expansion introduces errors in the strengths of all relevant, 
marginal, and irrelevant operators; all of which contain functions of longitudinal momenta 
because longitudinal locality is not maintained in light-front field theory Q. It has been 
shown that the hamiltonian computed to second order contains a two-body logarithmic 
confining interaction that acts in every Fock sector to confine quarks and gluons H]. This 
interaction alone is not rotationally invariant, and may require rotational symmetry to 
emerge dynamically only after different Fock sectors mix in the eigenstates. 

The largest errors in the final results may well originate from this perturbative approx- 
imation of the effective hamiltonian. That such an approximation may not be absurd is 
indicated by lattice gauge theory, where perturbative renormalization and tadpole resum- 
mation leads to extremely accurate results for heavy quark systems Q . Since errors in the 
relevant operators exponentiate in the scaling regime, it is likely that these operators must 
be fine-tuned. In lattice field theory one avoids most relevant operators by maintaining 
gauge invariance, although a relevant operator must be tuned to restore chiral symmetry 
if Wilson fermions are used. The relevant operators of immediate interest to us are quark 
and gluon self-energies, and in this paper we focus on the possibility that fine-tuning the 
gluon self-energy may be required to restore rotational invariance. 

The second step in our approach is solving for the hadronic bound states. If one tries 
to diagonalize the entire effective Hamiltonian, the wave function of a hadron must con- 
tain arbitrarily many parton components because the effective Hamiltonian still contains 
emission and absorption, albeit energetically limited. Instead, we divide the effective 
Hamiltonian into a part, Hq, which is solved nonperturbatively, and the remaining part, 
V, that is treated in bound state perturbation theory (BSPT). The division is arbitrary, 
but the choice of Hq determines the convergence of the calculation, and we must choose 
Hq to approximate the physics relevant for hadronic bound states. 
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If a constituent approximation emerges, we may hope to employ only interactions 
that preserve parton number in Hq. This is not required in our approach, since it is 
possible that mixing between few-parton sectors is non-perturbative without requiring 
many-parton sectors. If it is possible to completely ignore emission and absorption of 
partons in Hq, mesons can be approximated as quark-antiquark bound states with no glu- 
ons. However, rotational symmetry is maintained in this case only if the quark-antiquark 
potential is rotationally symmetric. Of course, we expect excited states in which the pair 
is separated by a large distance to contain gluons, so that it is only the short and inter- 
mediate range parts of the potential that need to be rotationally symmetric. Even in these 
ranges the potential need only be approximately invariant, since many-body effects found 
in higher orders of BSPT are required to restore complete rotational symmetry, as is well 
illustrated by light-front QED calculations [^. 

In previous work 1^, ^ a nonrelativistic reduction of the second-order effective Hamil- 
tonian was used to study mesons containing at least one heavy quark, approximated in 
leading order as color-singlet qq bound states. Hq included the short range Coulomb inter- 
action and a rotationally invariant long-range logarithmic potential obtained by angular 
averaging the long-range potential in the hamiltonian. V included a non-rotationally in- 
variant long-range potential, as well as low energy gluon emission and absorption, and 
other operators. The non-rotationally invariant potential leads to corrections in first-order 
BSPT, reported in |^ to be of order 30%. However, a mistake was found in the calculation 
of the non-rotationally invariant term. The calculation with the corrected potential re- 
veals that the corrections are actually less than 10%. Such corrections are not necessarily 
a serious problem, because second-order BSPT involving gluon emission and absorption 
will also lead to corrections of the same order when we use couplings required to produce 
heavy meson spectra. However, the calculation of these second-order BSPT corrections 
requires the solution of a three-body problem with confining two-body interactions, and 
we have not yet attempted this calculation. Since there are several ways in which rota- 
tional symmetry may emerge, we have chosen to investigate the simplest possiblity in this 
paper. 

The QCD hamiltonian contains quark and gluon self-energies that depend on longi- 
tudinal momenta even at second order. These self-energies and the logarithmic confining 
interaction produce a mass gap between the qq and qqg sectors, making a constituent 
approximation plausible. The self-energies are relevant operators, and the second-order 
approximation is probably inadequate; so we assign low- and intermediate-energy gluons 
an effective mass Mg{xg), and partially integrate them out. Note that we must allow 
this operator to be a function of the gluon longitudinal momentum fraction, Xg. This 
leads to new effective interactions in the Hamiltonian. The functional form of Mg{xg) is 
determined by the requirement that the violation of rotational symmetry in the effective 
interactions is reduced. We have found that this is achieved by using a x^-dependent 
effective gluon mass of the form Mg{xg) = rUgx'^. No other power of Xg, including a 
constant, or a polynomial in which the lowest term is not x^, leads to such an effect. 
By kinematical considerations, the parameter rUg is expected to be of order the glueball 
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massj^ In this letter we present numerical results obtained by applying the approach to 
heavy quarkonia for which the approximations are best justified. 

In previous applications |§, |3|, ^ the similarity transformation has been performed 
perturbatively around the free Hamiltonian. The starting point is a "bare" QCD hamil- 
tonian, H^, in the light-front gauge (^4+ = 0) regulated by a large cutoff , A, on energy 
differences at the interaction vertices. Counterterms must be added to this, -ff^'*'; and 
coupling-coherence fixes all couterterms. The transformed Hamiltonian Hx is expressed 
as 

Hx = SHaS^ = s[hI + hi') S\ (1) 

where 5 is a unitary transformation (5^5" = 1) such that the Hamiltonian Hx is band 
diagonal with respect to the scale A, i.e., 

Hx = fxGx , (2) 

where fx is a cutoff function which makes the Hamiltonian vanish outside the band. 
Equivalently, 

dHx 
d\ 

where 



[Hx,Tx] , (3) 



Ta = . (4) 

Let be the transverse cutoff (with dimension of mass^), a longitudinal scale (needed 
for dimensional reasons) , and let Eq denote eigenvalues of the free light-front Hamiltonian 
h. We can choose the cutoff function fx^^ = — |£'oi — E^j), for example, which 
fixes the transformation in Eq.(3). We then use Eq. (3) to determine how all operators in 
the Hamiltonian transform, and coupling coherence fixes the initial strengths of all non- 
canonical operators in terms of the canonical coupling, leading finally to an expansion for 
the renormalized Hamiltonian H^ in powers of the canonical coupling. 

A perturbative error is made in the strength of each operator, and this error typically 
grows as A decreases because the canonical coupling grows. As stated above, errors in 
the strength of relevant operators are expected to produce the largest bound state errors, 
and we assume that at some scale these errors become sufficiently large that we must fix 
them by adding mass terms that are tuned by hand. We will add a gluon mass term at a 
scale A, which from this point will refer only to the specific cutoff where the gluon mass 
is added. Moreover, since the gluon mass makes a significant contribution to the diagonal 
components of the Hamiltonian, we alter the transformation so that we now perturb about 
the free Hamiltonian with the gluon mass included. The first step is to lower the cutoff 
to the scale A and obtain the renormalized Hamiltonian, 

H^ = h + fxvx- (5) 



^Note that nig is not one-half the glueball mass because of the a:-dependence in this operator. 
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At this point we add an effective a;-dependent gluon mass operator to and regroup 
the Hamiltonian so that the new one-body operator is included in the free Hamiltonian 

H^^H'^ = h' + hvx. (6) 

i •'I 1 

Below the scale A the new dispersion relation for gluons is given by: 

^ -It — — • 

Now, we run a second similarity transformation. Since we have changed the free Hamil- 
tonian, we have to use a similarity transformation that runs a cutoff for the new gluon 
dispersion relation. This second transformation is implemented by introducing a new 
similarity function, f\. 

At this point the Hamiltonian is band diagonal with respect to the cutoff function /, 
but not with respect to the cutoff function /. We cannot simply change the cutoffs by 
hand, so we must run the second transformation from infinity to whatever final value of 
the cutoff we choose. The procedure for running this second transformation is formally the 
same as that for running the first transformation. There is no need for new countcrtcrms, 
since we have already completely renormalized the hamiltonian and the second transfor- 
mation does not generate any new divergences. The output of the second transformation 
is an 0{g^) effective Hamiltonian band-diagonal in light-front energy with respect to the 

2 

scale and can be written as : 

= + ^(1) + ^(2) + + . (9) 

The light-front kinetic energy is augmented by the effective gluon mass operator, and 
it is still diagonal. The 0{g) term v^^^ now corresponds to (massive) gluon emission 
and absorption with nonzero matrix elements only between states with energy difference 

2 

smaller than The gluon mass sets a dynamical threshold for gluon emission and 

absorption. Since the mass of a free gluon can be arbitrarily small depending on the Xg- 
dependence chosen, particle number changing interactions are, in general, not integrated 
out completely. This means that there still may be a contribution from the qqg bound 
state, but compared to the previous work where the low-energy gluons were left massless 
in the effective Hamiltonian, this contribution is significantly reduced. The 0{g'^) canon- 
ical instantaneous gluon exchange interaction, v^'^\ and the 0{g^) effective interactions 

(2) 

generated by the first transformation, vlff, remain unchanged, apart from the new overall 
similarity function that makes the Hamiltonian band diagonal with respect to the scale 
The new 0{g'^) effective interactions, v^^ff^ generated by the second transformation 
depend on the gluon mass rUg and its functional form. 
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Using Mg{xg) = mgx'^ the effective qq potential in the nonrelativistic limit, including 



the new terms is given by: 



^-\Eoi-Eoj\Ug^CFi2Mf 




A2 q2(2M)^ 



+ 



gi 



1z q2 + 



(2M)2 



' q2(2M)\ 
,7'+ \qz\P+ I 



V 



\qz\P+ 




^ > (10) 



where q = {qz,q±) is the exchanged (gluon) momentum, M is the quark mass, Cp = 
4/3 is a color factor and is the meson total longitudinal momentum. To take the 
nonrelativistic limit, we introduce a z-component for momenta. 



pi + + M2 + 



(11) 



so that qz = piz — P2z, and we make an expansion in powers of p/M, dropping all non- 
leading terms. 

The potential in eqn. (10) results from instantaneous gluon exchange in the light- 
cone gauge added to terms from the two similarity transformations. Each similarity 
transformation produces a second-order gluon exchange interaction that replaces direct 
exchange of gluons not allowed by the cutoffs on energy change at the vertices. If we had 
not added a gluon mass and run a second transformation, the parameters would be the 
cutoff /i, the quark mass M, and the coupling a. The second similarity transformation 
introduces two additional parameters, A and mg. 

The purpose of this work is to explore the possibility that the gluon mass operator may 
lead to an effective quark-antiquark potential that is more nearly rotationally invariant 
than the potential computed in our previous work. To show this, we choose Hq as in 
our previous work (i.e. the nonrelativistic reduction of the kinetic energy, the effective 
one-body operators, Coulomb potential and the angular average of the confining potential 
with constituent masses). The bound state calculation proceeds in full analogy with [Q, ^. 
We assume that the scale n is small enough that the lowest-lying states are dominated by 
their qq component. 



The potential given in Eq. |10| is Fourier transformed and expanded in even Legendre 
polynomials, 

oo 

V{r) = Ko,i(r) + ^ yconf,2k(r-)P2k(2) , 

k=0 

where r is the quark-antiquark separation, z = cos{9) and 



(12) 



V'conf,2k(?') 



dz y(r)P2k(^) 



The eigenvalue problem for the leading order hamiltonian Hq is 

Ho\P)=M^\P) , 



(13) 



(14) 
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where Ai'^ is the invariant mass of the bound state and A4'^ 
This leads to a dimensionless Schrodinger equation 



4M2 + 4ME defines E. 



-4r + C (-Vconifii^) + ^coul(7^)') 



with 



7^ = £r, 



2maCp 



2m{E - S) 
£2 : 



(15) 



(16) 



where C = {fx^ /V^){P^ /2M) , m = M/2 is the reduced mass, a = g'^/'iir and S is 
the finite shift produced by the self-energies after subtracting terms needed to make the 
confining potential vanish at the origin, 



TT 
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A2 



In 



4M2 



ml + 4M2 



(y2 + u;) ^ m2y2 



y(l-y) 4M2y 



1-y 



y2 + tt; + y2(-x _ y)m2/(4j^^2) ' 



(17) 



The parameters are adjusted so that the angular average of the potential is close to the 
one used for charmonium in our previous work, as displayed in Fig. 1. For a qualitative 
analysis this is sufficient. Here, the potential obtained in ref. is calculated by using 
M = 1.6 GeV, ji = 1.8 GeV and a = 0.53f| For the new potential we use M = 1.6 GeV, 
rUg = 1.6 GeV, A = 9.5 GeV, ^ = 1.5 GeV and a = 0.55 to obtain a reasonable fit for the 
masses of the lowest lying states of charmonium: A^i^ = 3.1 GeV, Mip = 3.5 GeV and 
M2S = 3.6 GeV. 

In Fig. 2 we compare the angular dependence of the full potential to the one obtained 
in our previous calculation. The violation of rotational symmetry in the new effective 
potential is significantly reduced in the intermediate range (0.3 — 1.5 fm) and increases 
at large distances. As stated above, the long range part of the potential need not be 
rotationally invariant because states in which the qq pair are separated by large distances 
will contain valence gluons. At short distances (< 0.3 fm), where the Coulomb potential 
dominates, the violation is very small, although larger for the new potential. With the 
parameters fixed we evaluate the corrections due to rotationally noninvariant terms in 



^An error in the expression for the self energy in ref. |^ has been corrected, 
parameters required to fit the lowest lying states slightly from those in . 



This changes the set of 
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the confining potential. There are corrections to P-states in first-order BSPT, while S- 
states are corrected in second-order. For weak couplings these corrections scale with the 
coupling at the same order as second-order BSPT corrections involving gluon emission 
and absorption. We concentrate on the corrections to the lowest lying P-state, which are 
easily computed. These corrections to the lowest lying P state of heavy quarkonium are 
reduced from about 8% |^] to about 2%. 

In conclusion, we have shown that it is possible to adjust the relevant gluon mass op- 
erator so that apparent violations of rotational invariance in heavy quarkonia are reduced. 
It is assumed that non-perturbative corrections to this operator become important as the 
hadronic bound state scale is approached by the cutoff, and that fine-tuning this operator 
allows us to perturbatively reduce the cutoff further without modifying other operators by 
hand. We use rotational invariance as a test of these ideas, tuning the dependence of the 
gluon mass on longitudinal momentum to reduce naive violations of rotational invariance 
by a factor of 3 — 4. Rotational invariance of the two-body interaction is significantly 
improved at intermediate distances if a gluon mass operator of the form rn^x^ is used. 
Corrections to rotational invariance are sufficiently small in charmonium even without 
such a gluon mass that we can not yet conclude that this operator is required. A real 
test of the significance of the relevant gluon mass operator requires higher order BSPT 
calculations and a study of light mesons. 
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References 

[1] K. G. Wilson, T. S. Walhout, A. Harindranath, Wei-Min Zhang, R. J. Perry and 
S.D. Glazek, Phys. Rev. D 49, 6720 (1994), |hep-th 7940115^ . 



[2] R. J. Perry, in Proceedings of Hadrons 94, edited by V. Herscovitz and C. Vasconcellos 
(World Scientific, Singapore, 1995), |hep-th/9407056t R. J. Perry, in Proc. Fourth Int. 



Workshop on Light-Front Quantization and Non-Perturbative Dynamics, edited by 
S. D. Glazek, (World Scientific, Singapore, 1995), p.56, |hep-th/9411037t 



[3] M. Brisudova, R. Perry and K. Wilson, Phys. Rev. Lett. 78, 1227 (1997). 
[4] K. G. Wilson, Rev. Mod. Phys. 55, 583 (1983). 

[5] St. D. Glazek and K. G. Wilson, Phys. Rev. D 48, 5863 (1993); ibid. 49, 4214 (1994); 
F. Wegner, Ann. Physik 3, 77 (1994). 

[6] R. J. Perry and K. G. Wilson, Nuc. P hys. B 403, 587 (1993); R. J. Perry, Ann. 
Phys. 232, 116 (1994), |hep-th/94020T5i Similar ideas were previously developed by 



R. Oehme, K. Sibold and W. Zimmerman, Phys. Lett. B 147, 115 (1984). 



8 



[7] 



[8] 
[9] 



G. P. Lepage, Talk given at 35th Internationale Universitatswochen fuer Kern- und 
Teilchenphysik (International University School of Nuclear and Particle Physics): 
Perturbative and Nonperturbative Aspects of Quantum Field Theory, Schladming, 
Austria, 2-9 Mar 1996. In Schladming 1996, Perturbative and nonperturbative aspects 
of quantum field theory 1-48, |hep-lat / 9607076 ; 

Lectures given at 9th Jorge Andre Swieca Summer School: Particles and Fields, Sao 
Paulo, Brazil, 16-28 Feb 1997, |nucl-th/9706029 . 

M. Brisudova and R. Perry, Phys. Rev. D 54, 1831 (1996), |hep-ph/95H"443 . 



B. Jones, R. Perry, S. Glazek, Phys. Rev. D 55, 6561 (1997), hep-th/9605231 



9 



Figure Captions 



Figure 1: The angular average of confining and Coulomb potential witli massive gluons 
compared to the angular average with massless gluons. The parameters of the potential 
with massive gluons were adjusted so that it is close to the potential with massless gluons. 

Figure 2: The angular dependence of the potential with massive gluons compared to the 
potential with massless gluons at r = 0.2 fm, 1.0 fm, 2.0 fm. The violation of rotational 
symmetry is significantly reduced in an intermediate range (0.3 — 1.5 fm). At very short 
distances the violation of rotational symmetry is small, although larger for the potential 
with massive gluons. 
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Figure 1: The angular average of confining and Coulomb potential with massive gluons com- 
pared to the angular average with massless gluons. The parameters of the potential with 
massive gluons were adjusted so that it is close to the potential with massless gluons. 
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Figure 2: The angular dependence of the potential with massive gluons compared to the poten- 
tial with massless gluons at r = 0.2 fm, 1.0 fm, 2.0 fm. The violation of rotational symmetry 
is significantly reduced in an intermediate range (0.3 — 1.5 fm). At very short distances the 
violation of rotational symmetry is small, although larger for the potential with massive gluons. 
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